Non-Markovian entanglement dynamics in the presence of system-bath coherence 
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A complete treatment of the entanglement of two-level systems, which evolves through the contact 
with a thermal bath, must include the fact that the system and the bath are not fully separable. 
Therefore, quantum coherent superpositions of system and bath states, which are almost never 
fully included in theoretical models, are invariably present when an entangled state is prepared 
experimentally. We show their importance for the time evolution of the entanglement of two qubits 
coupled to independent baths. In addition, our treatment is able to handle slow and low-temperature 
thermal baths. 



The crucial feature of quantum information is the pres- 
ence of coherent superpositions. A single isolated two- 
level system (qubit) can be prepared in a superposition 
of its and 1 states, and the manipulation of such states 
leads to new possibilities for the storage and processing 
of information. In a pure quantum system, the superpo- 
sition is entirely coherent, which means that a definite 
phase relation exists between the and 1 states. Unlike 
in the ideal isolated case, the interactions of real quantum 
systems with their environment lead to the destruction 
of these phase relations, in other words, decoherence. 

More interesting than a single two-level system is a 
collection of multiple such qubits. Coherent superposi- 
tion states of multiple qubits can be prepared, and their 
dynamics are studied to understand the functioning of 
quantum networks X The presence of phase relations be- 
tween qubits, a second type of coherence, is termed en- 
tanglement. The destruction of entanglement through 
interactions with the environment is an important prob- 
lem, both for a fundamental understanding of quantum 
mechanics and the development of quantum information 
processing^ Recently, it has been found that the decay of 
entanglement can be very different from that of the single 
qubit coherence. In particular, the entanglement can dis- 
appear completely in finite time^ Interactions between 
qubits further affect the time scale of the decay^ 

Here, we argue that, besides single qubit coherence 
and the entanglement of multiple qubits, there is a third 
form of coherence that is important for the description 
of quantum information. The physical nature of this ef- 
fect lies in a detailed description of the heat bath that 
leads to decoherence and disentanglement. In a complete 
theory, both the system and the bath are quantum me- 
chanical. It is therefore possible, and often unavoidable 
in experiment, to create coherent superpositions of sys- 
tem and bath states, and such superpositions will play an 
important role in the time evolution of the qubit system 
of interest. 

The description of entanglement dynamics starts from 
an equation of motion for the qubits together with an 
initial condition, which are both affected by a correct 
description of the bath. Most theoretical treatments, 
including the Rcdficld and Lindblad formulations, de- 
scribe the qubit dynamics under the Born and ultrafast 
bath approximations. 5 These lead to convenient time- 



local equations of motion, but are only valid if the time 
scale of bath dynamics is much faster than the charac- 
teristic time scales of the system.— Furthermore, in the 
same limit, the initial state can be chosen to be of a form 
in which system and bath are independent. If the fast 
bath condition is not fulfilled, however, the approxima- 
tion breaks down. As briefly explained below, it is also 
invalid at low temperature, where quantum fluctuations 
contribute longer time scales^ 

Although the effect of noise correlation (or non- 
Markovian dynamics) on the entanglement has been 
studied under the rotating wave approximation at zero 
temperature^ the role of the bath on the initial condi- 
tion has rarely been investigated. However, in the regime 
where non-Markovian effects are important, the presence 
of system-bath correlations invalidates the initial state 
in which the system and the bath are independent. Es- 
pecially in the experimentally relevant case where the 
qubit system is excited out of equilibrium and the sub- 
sequent dynamics is probed, a proper treatment of the 
initial state is crucial. 

It is the role of the bath in the equation of motion and, 
in particular, in the initial state, that we explore in this 
paper. We present a new method to calculate the entan- 
glement that rigorously deals with the quantum dynam- 
ics in this situation. The method treats the system-bath 
interaction non-perturbatively and without assuming a 
fast noise bathiiir— It thereby enables us to study the 
role of system-bath coherence and of an initial state in 
which a coherent superposition of the qubit system and 
the bath is prepared. 

Before introducing our rigorous approach, we clarify 
the role of system-bath coherence by starting from a 
simple model, in which both the system and the bath 
are two-level systems with transition energy e. In this 
much simplified situation, the time scales of system and 
bath evolution are identical, and therefore the role of 
system-bath coherence is expected to be strong. With 
a coupling g the complete Hamiltonian is given by H — 
ecgCs + tc^cs + <?(cgCB + c B cs), where c$/b and cL B are 
the usual annihilation and creation operators working on 
the system or bath. The system plus bath are assumed 
to be in thermal equilibrium at an inverse temperature 
j3. Apart from a shift in the effective excitation energy, 
no effect of the bath is apparent in the reduced density 
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matrix for the system degrees of freedom p. By consider- 
ing the full density matrix R of system plus bath, how- 
ever, one finds a matrix element that represents coher- 
ence between the system and the bath, (Ogl b|-R|1s0b) = 
— ^e _/3c smh/3g, where S denotes the system and B the 
bath, and Z is the total partition function. This term 
would be zero for the commonly considered initial con- 
dition R = p® exp(— PHb)/^ exp(— /3Hb)- Its presence 
shows that the density matrix in equilibrium is not sep- 
arable; it cannot be written as the product of a system 
and a bath part. The correct treatment of the system- 
bath coherence in the initial state is important for the 
entanglement dynamics following a pulse that excites the 
system out of thermal equilibrium. 

In the following, we will study the entanglement of two 
qubits each separately coupled to a more general bath, 
which evolves on a characteristic time scale 1/7. The 
qubits are labeled 1 and 2, both have an excitation energy 
e, and they are coupled by an interaction J. The system 
Hamiltonian is 



e(c\ci + c\c 2 ) + J{c\ + ci)(4 + c 2 ). 



(1) 



To study system-bath coherence, we need to introduce 
a fully quantum-mechanical bath. A bath that is suffi- 
ciently general to model many physical systems, while at 
the same time allowing efficient calculations, is given by a 
set of harmonic oscillators^ The individual bath modes 
are labeled with an index j, and have masses rrij, fre- 
quencies ujj , coordinates Xj and corresponding momenta 
Pj . The Hamiltonian for the system-bath interaction plus 
the harmonic bath is (fi. = 1 throughout the paper) 



HsB+B — — y^ffqjVq" 




(2) 



Starting from this Hamiltonian, we need to derive an 
equation of motion for the reduced density matrix of the 
system. It must describe the fluctuations in the system 
energies and the dissipation arising from the interaction 
with the bath. Furthermore, we want the method to be 
able to describe an initial state in which the system and 
bath are correlated, i.e. for which the complete density 
matrix cannot be written as a direct product of a system 
and a bath part. 

In the reduced description, all necessary information 
on the bath and the system-bath coupling is contained 
in the spectral densities J a(x >(uj) — J2j ^"li' 3 ~ ^i )- 
While the V Q 's can be any system operator, we focus 
here on two independent baths which induce transitions 
in a single qubit each, given by V a = c a +cj, for a = 1, 2. 
Thus, the spectral representations of the cross-correlation 
functions are identically zero, J\2{u) = 3i\(y>} = 0. Fur- 
thermore, we assume that each bath evolves on the same 
time scale and couples to a qubit with the same strength. 
Both baths can then be described with the same spectral 
density function, Ju{uj) = ^22^) ='■ Ji^), which we 
model as J(u) = uj^jj. To see the effect of fluctu- 
ations in the system parameters and of the dissipation 



originating from this spectral density, we consider the 
correlation function of the noise variables. It can be cal- 
culated as the Fourier transform of the spectral density 
to be 
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where we have defined the time scales 1/vk — /?/27rfc 
(for k > 1) and 1/vo = 1/7 and the prefactors Co = 
A7H + cot/37/2), c fc = (4A7//3)i/ fc /(i/^ - 7 2 ). For large 
7, that is, in the case of a fast bath, the dissipative term, 
given by the imaginary part of the correlation function, 
becomes delta correlated. Quantum effects, however, 
make the fluctuations induced in the system (given by 
the real part of L{t)) always correlated for low tempera- 
ture, in which case the evolution of the bath contains time 
scales fixed by the Matsubara frequencies i^.— The fast 
bath (Markovian) approximation is then invalid. Note 
that we do not make the rotating wave approximation in 
the system-bath interaction, which can lead to significant 
changes in entanglement dynamics^ 

Using the Feynman- Vernon influence functional ap- 
proach, the reduced equation of motion can be written in 
the form of path integrals over the system coordinates. 
The path integrals can be evaluated numerically ) 14 i 15 or 
the equation of motion can be rewritten as a stochastic 
Schrodinger equation.»A£ Alternatively, a time-local equa- 
tion of motion for the reduced density matrix can be 
derived^— It accounts for the dynamics of the quan- 
tum bath and the system-bath coherence through a set 
of auxiliary density matrices, whose time evolution is also 
governed by time-local equations of motion. The set of 
density matrices can be propagated efficiently in a com- 
puter, and this description is therefore well-suited for the 
study of entanglement dynamics. The complete equation 
of motion involves multiple auxiliary density matrices, la- 
beled by a multi-index n Q fe, where a labels the various 
terms in the system-bath interaction, and k the Matsub- 
ara frequencies included to allow for low temperature. It 
includes the time scales and prefactors from the correla- 
tion function, and is given hyl*^ 
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where A x B = [A, B]. The notation n a k — > n a k + 1 
refers to an increase in the afc'th component of the multi- 
index, while all other indices are unchanged. Similarly, 
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FIG. 1: Decay of entanglement from an initially maximally 
entangled state \<f>) = (|1) — |2))/v2. The solid line is the 
result from our full calculation while the Redfield result is 
plotted as a dash-dotted line for comparison. 

p n a k = p n ak^n ak -i denotes a decrease of this index. For 
the simulations in this paper, we will use the values of e = 
1.5J, A = 0.3J, 7 = 0.5J, (3 = 2.5J. This method has 
recently been applied to calculate the exciton dynamics 
in light-harvesting complexes.— ~— 

Given the equation of motion, the dynamics of the sys- 
tem can be calculated numerically starting from an arbi- 
trary initial condition. The procedure takes the system- 
bath coupling into account without invoking perturba- 
tive, Markovian, or rotating wave approximations. The 
system-bath coherence in the initial state can be included 
through the auxiliary density matrices. In particular, to 
obtain a thermal equilibrium initial state, which will in 
general be a partially coherent superposition of the sys- 
tem and the bath, one can propagate the equation of mo- 
tion starting from any state for a time longer than any 
characteristic time scale in the system or the bath. Here, 
we will not consider the possibility of a non-ergodic sys- 
tem, in which case additional averaging is required. We 
will now first discuss the effect of the bath in the equa- 
tion of motion, followed by a treatment of the correlated 
initial state. 

A popular measure of the entanglement of two qubits 
is the entanglement of formation, which can be calcu- 
lated directly from the reduced density matrix through 
the concurrence 3i The effect of slow bath dynamics has 
been shown to fundamentally alter the time evolution of 
entanglement between two quantum systems at zero tem- 
perature. In particular, strong coupling to a bath can not 
only change the time scale on which entanglement disap- 
pears, but also lead to revival of entanglement after a 
period of zero concurrence^ This effect is also observed 
in the current model. In FigfTJ the concurrence is plot- 
ted as a function of time, starting from an initially en- 
tangled state |0) = (|1) — \2))/y/2, and the ad-hoc but 
common assumption of an independent bath. The dash- 
dotted line shows the result of a Redfield calculation in 



the secular approximation, which treats the system-bath 
interaction perturbatively and in the Born and Markov 
approximations^^ Without the secular approximation, 
the Redfield equation can lead to an unphysical den- 
sity matrix, and the concurrence is ill-defined. The re- 
sult from the full calculation (shown as a solid line) is 
markedly different from the Redfield predictions. Before 
reaching the same positive equilibrium concurrence as in 
the Redfield calculation, it shows a revival of concurrence 
(at t f=a 3) after sudden death of entanglement (around 
t sa 2). This calculation shows that conventional quan- 
tum master equations do not completely describe the en- 
tanglement dynamics. To understand the positive equi- 
librium value of the concurrence, which is the combined 
effect of the coupling between the two qubits and the 
temperature, we consider the following simplified model. 

For two qubits weakly interacting with a heat bath, one 
would expect the reduced density matrix in equilibrium 
to be given by the Gibbs measure p cq = e~^ Hs /Tre~^ Hs . 
For a system Hamiltonian Hg — e(c{ci + c^) + J(c\c2 + 
cjci), the concurrence in this equilibrium state is 

= -1 + sinh/3J 

cosh/3e + cosh/3J' 1 ' 

which is positive for j3J > sinh~ 1 ss 0.88. We thus find 
that for large J or at low temperatures, the equilibrium 
concurrence is positive. This suggests that it might be 
better to consider the entanglement of the eigenstates 
of the system Hamiltonian instead, which vanishes in an 
equilibrium state obtained through weak coupling with a 
heat bath. 

This is only half of the story. Apart from its effect 
on the equation of motion, the slow bath also makes it 
impossible to factorize the initial state into a system and 
a bath part. This becomes clear if we consider a realis- 
tic physical initial state. In the previous discussion, the 
initial state was assumed to be a product state of system 
and bath density matrices. This is non-physical in the 
case of strong system-bath coupling, because it is usu- 
ally hard, if not impossible, to prepare such a state. 

A more realistic initial state can be obtained by ap- 
plying a pulse to the equilibrium state. Here, we con- 
sider a pulse that affects only the first qubit by rotating 
it over an angle n. This operation can be expressed as 
p — » a\ pa\* . Such an operation clearly takes the system 
out of equilibrium, and it creates additional entangle- 
ment. It can be realized for example as a 7r-pulse in nu- 
clear magnetic resonance^ The resulting time evolution 
of the concurrence is shown in Fig. [5] as a solid line. The 
entanglement dynamics is not captured in the Redfield 
approach (which cannot include the system bath corre- 
lations in the initial state), shown as a dash-dotted line. 
To see the effect of system-bath coherence explicitly, we 
also perform the calculation for a factorized initial state 
(resulting in the dynamics plotted as a dashed line in 
Fig. [5]). Although the reduced density matrices for both 
the system and the bath are the same as in the full cal- 
culation, the system-bath coherence was set to zero just 
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FIG. 2: Time evolution of the entanglement after a pulse has 
been applied to the equilibrium state. The solid line shows the 
result of the full calculation, while the system-bath coherence 
at t = was neglected in the results presented as a dashed 
line. The dash-dotted line shows the Redfield result. 

before applying the spin flip operation at t = 0. Even 
though the system-bath coherence is fully accounted for 
in the subsequent time evolution in both cases, it is clear 
that the presence of the system-bath coherence in the 
initial state changes the entanglement dynamics dramat- 
ically. The initially created entanglement is much larger 
than in the case of a product initial state. Moreover, the 



effect lasts long enough to reduce the entanglement death 
time considerably. Death of entanglement occurs only for 
short periods of time in the full calculation, while no en- 
tanglement is found until after t = 4 when the initial 
system-bath correlation is ignored. This finding shows 
that, even in calculating a property such as the entan- 
glement that depends only on the system degrees of free- 
dom, the initial system-bath coherence must be properly 
taken into account. In effect, memory of entanglement 
can be stored in the bath. While we have uncorrelated 
baths coupled to each qubit in this study, the effect of 
correlation in bath modes might well further strengthen 
this finding, and entanglement measures that include the 
system-bath coherence will be needed. 

In conclusion, we have studied the entanglement of two 
qubits in the presence of a quantum mechanical bath. 
The coherence between system and bath is found to have 
an important effect on the time evolution of the con- 
currence when the system is excited out of equilibrium, 
which is not correctly described by a conventional quan- 
tum master equation. These findings will be relevant 
for the design of quantum networks 1 and information de- 
vices, as well as for the dynamics of excitations in bio- 
logical systems ) 17 ' 18 ' 23 quantum dots^i and conjugated 
polymers^ For spin systems, the effect of a fermionic 
bath should be investigated. 26 

We thank dr. A Ishizaki for helpful comments on the 
manuscript. 
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